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I. INTRODUCTION 



It is theorized that Einstein's General Relativity breaks down at sufficiently high energies. 
Therefore, General Relativity may not be the most appropriate choice for a gravitational 
theory in the very early universe. Developments in string theory suggest that gravity may 
truly be a higher dimensional theory, becoming an effective 4-dimensional theory at lower 
energies. Some researchers have suggested an alternative scenario in which the matter fields 
are restricted to a 3-dimensional brane- world embedded in 1 + 3 + d dimensions (the bulk), 
while the gravitational field is free to propagate in the d extra dimensions [0, |^. This 
new scenario for the gravitational field will have significant implications for cosmology. 

The dynamical equations on the 4-dimensional brane- world differ from the field equations 
of General Relativity, since there exist additional terms that carry non-local gravitational 
effects from the bulk onto the brane. These additional terms are quadratic in nature and 
reduce to the regular Einstein Field equations of General Relativity for late times. How- 
ever, for early times, these quadratic terms will play a very critical role in the evolutionary 
dynamics of these brane-world models. 

A lot of effort has been directed at the so-called perfect fluid Friedmann Brane-world 
models P, H, ^ H, 0, |[. More recently, people have begun to investigate anisotropic brane- 
world models |^, ||, |^, |T0|, |Tl], |T2|, trying to determine the effects of the bulk gravitational field. 
People have also begun to look at the effect of scalar fields in the Friedmann brane-world 
p!3| . It is known in General Relativity that only sets of measure zero in the set of all spatially 
homogeneous models (satisfying a set of reasonable energy conditions) will isotropize []T^ 



It is not yet clear whether this result holds for spatially homogeneous brane-world models 
as well. 

General relativistic cosmological models containing both a perfect fluid and a scalar field 
having an exponential potential have been previously studied |16|, One of the exact 
solutions found in these models is a spatially flat isotropic model with the property that the 
energy density of the scalar field is proportional to the energy density of the perfect fluid, 
and consequently are labeled scaling cosmologies (See |]T8|, |T9[ and references therein). It 



was found in |jT5|, |T6|, that these cosmological scaling solutions are in general not stable, 
and hence do not represent typical late-time behaviour. However, during the evolution of 
one of these cosmological models containing both a perfect fluid and a scalar field having an 
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exponential potential, these scaling solutions may play an important role in the intermediate 
or transient behaviour of these cosmologies, since the scalar field and the perfect fiuid are 
non-negligible in a neighborhood of the scaling solution. In particular, in these models a 
significant portion of the current energy density of the Universe could be contained in the 
scalar field whose dynamical effects mimic those of cold dark matter. 

In general relativistic cosmological models, the Bianchi II model plays an important 
role in the past asymptotic behaviour of the Bianchi IX models. Recently, Coley 0] has 
argued that the Bianchi type IX model has an isotropic singularity in the brane-world 
scenario. This indicates that the Bianchi II model in the brane-world scenario may also have 
a behaviour that is different than in the general relativistic case. In this paper, we shall 
assume that the brane-world is a spatially homogeneous Bianchi type II model containing 
both a scalar field and a perfect fiuid. Consequently, due to the additional terms that carry 
non-local gravitational effects from the bulk onto the brane, new evolutionary behaviours 
in the cosmologies are expected. In particular, with a combination of a perfect fluid and a 
scalar field, due to the quadratic nature of the non-local gravitational effects, a wide variety 
of new and exotic asymptotic and transient behaviours are possible. 

Maartens and Shiromizo et al. [pT| , |22| have developed an elegant covariant approach 
to the bulk effects on the brane. The equations derived by Maartens are an extension of 
earlier work by Ellis and MacCallum that has been subsequently developed more recently 
in the book by Wainwright and Ellis Using the formalism developed by Maartens and 
Wainwright and Ellis, we propose to investigate the dynamical behavior in a wider class of 
anisotropic models than what has been previously analyzed. 

The resulting field equations will yield a system of ordinary differential equations, suitable 
for a geometric analysis using dynamical systems techniques. This analysis will determine 
whether the dynamics of the brane-world scenario mimics the dynamics of a General Rela- 
tivistic cosmology or at the very least, is asymptotic to a general relativistic cosmology. One 
would be interested in both the early (nature of the initial singularity) and late time behavior 
(i.e., whether these models approach the General Relativistic regime and isotropize). 

With regards to notation, quantities having a tilde are objects found in the five dimen- 
sional bulk while non-tilde quantities are objects restricted to the four dimensional brane- 
world. The upper case latin subscipts {A, B,C, . . .) range from 0..4, the lower case latin 
subscripts a, 6, c, . . . ) range from 0..3 and the greek subscripts (a, /5, 7, . . . ) range from 1..3. 



3 



II. GOVERNING EQUATIONS 



The popular brane-world scenario assumes the existence of a domain wall, M, (a 3-brane) 
with induced metric qab = Qab — nAUB in a five dimensional spacetime M with metric gAB 
where is the spacelike unit vector normal to M. In a neighborhood of the brane, one can 
choose coordinates = {x'^, x) such that ua = X,a where x = coincides with the brane. 
The are the spacetime coordinates on the brane M. The five dimensional field equations 
are Einstein's equations with a bulk cosmological constant A, and a brane energy-momentum 
source Tab, 

Gab = K'^[-AgAB + S{x){>^gAB + Tab}] 

where k"^ = Sn/Mp, with Mp being the five- dimensional Plank mass and A is the brane 
tension. 

The field equations induced on the brane have been derived using a geometric approach by 
Shiromizu et al. ||2D|, using the Gauss-Codazzi equations, Israel's junction condition 

and imposing a Z2 symmetry with the brane as a fixed point. The result is a modification 
of the standard Einstein equations with new terms {Sab and Sab which are defined later) 
carrying bulk effects onto the brane: 

Gab = -Agab + K^Tab + 'f^^Sab " £ab (2.1) 

where 



= ^ A = 6— A 



A+|4^1A^ 

3Mp3 



(2.2) 



Ml k4 

with Mp being the four-dimensional Plank mass, and A is the effective four- dimensional 
cosmological constant on the Brane. It is common to assume that the effective cosmological 
constant. A, on the brane is zero through fine tuning of the brane tension A. However, we 
shall assume for the moment that it is non-zero. 

We will assume that the source term. Tab (restricted to the Brane) is a non-interacting 
mixture of ordinary matter having energy density p, and a minimally coupled scalar field, 
0, that is 

Ti rp perfect fluid . rp scalar field /n o\ 

ab-J-ab +-^ab ■ l^-^J 

We will assume that the matter component is equivalent to a co-moving perfect fiuid with 4- 
velocity u"" and a linear barotropic equation of state p = (7 — l)p. Energy conditions impose 



4 



the restriction p > and causality requires that 7 G [0, 2]. We will assume that the potential 
for the scalar field has an exponential form, that is, V{(f)) = Voe'"^'^ = Vqc ^h^'^. The energy- 
momentum tensors restricted to the brane for a perfect fiuid and a minimally-coupled scalar 
field are 

y^^pe./ect/W ^ pu^u,+pK,, (2.4) 
j^^^scalar freM ^ ^^^^^^ _ ^ ^ (2.5) 

where the projection tensor hah = dab + UaUb projects orthogonal to m". We note that if 0'" 
is timelike, then a scalar field with potential V{(f)) is equivalent to a perfect fiuid having an 
energy density and isotropic pressure 

^scalar field ^ _ 1^. c ^ y ^ (2.6) 
^scalar field ^ -l^.^^i^ - V {(f)) . (2.7) 

The bulk corrections to the Einstein equations on the brane are of two forms: firstly, the 
matter fields contribute local quadratic energy-momentum corrections via the tensor Sab, 
and secondly, there are nonlocal effects from the free gravitational field in the bulk. The 
local matter corrections are given by 

Sab = ^T^^'Tab - \TacT\ + ^^Qab [^T^dT^' - (T/)^] . (2.8) 

Since the Sab is essentially quadratic in the energy momentum tensor, we should expect 
cross terms between the scalar field and the perfect fiuid. In this paper we shall assume 
that the gradient of the scalar field 0'", is aligned with the fiuid 4- velocity, u"", that is 
0'"/ a/— = u""- In general 0'*^ need not be aligned with u"" thereby creating a rich 
variety of cross terms. The local brane effects due to a combination of a perfect fluid and a 
scalar fleld are 

1/1 



Sah = -^yp- ^(l>;c(i)''' + V{(f)) j UaUb 

+ ^ (p - ^0;c0^^ + V^(^p + 2p- ^0^,0^^ - y(0)^ hat (2.9) 

The non-local effects from the free gravitational fleld in the bulk are characterized by 
the projection of the bulk Weyl tensor onto the brane. Given a timelike congruence on the 
brane, the bulk correction. Sab can be decomposed pO[ via 



ab 



1 

U{UaUb + -hab) + Vab + 2Q(aMf,) 



(2.10) 



(See pO[ for further details.) In general, the conservation equations (the contracted Bianchi 
identities on the brane) do not determine all of the independent components of Sab- In 
particular, there is no equation to determine Vab and hence, in general, the projection of the 
5-dimensional field equations onto the brane do not lead to a closed system of equations. 
In the cosmological context, in which the background metric is spatially homogeneous and 
isotropic, we have that 

DaU=Qa=Vab = (2.11) 

where is the totally projected part of the brane covariant derivative. Since Vab = 0, in 



this case the evolution of £ab is fully determined |2^, |2g, ^ 



All of the brane source terms and bulk corrections mentioned above may be consolidated 
into an effective total energy density, pressure, energy flux, and anisotropic pressure as 
follows. The modified Einstein equations take the standard Einstein form with a redefined 
energy-momentum tensor: 

Gab = n'Ta,'°'^' (2.12) 

where 

(2.13) 



A 1 

rp total _ nJ-Tj- Q— P 

-L ab — riUab ' ab ~r r, *Jafe r,(^ab- 



The total equivalent energy density, pressure, energy flux, and anisotropic pressure due to 
a perfect fluid, scalar field, and both local and non-local brane effects are 
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(2.14) 
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(2.15) 

(2.16) 
(2.17) 



As a consequence of the form of the bulk energy-momentum tensor and the Z2 symmetry. 



it follows [^, |2^ that the brane energy-momentum tensor separately satisfies the conser- 
vation equations, (where we have tacitly assumed that the scalar field and the matter are 



non-interacting) that is, 



rpa perfect fluid ^ /o i o\ 



■ 6; a, 

Tia scalar field ^ /r) i n\ 

Given that T^"^ = 0, the contracted Bianchi identities on the brane imply that the projected 
Weyl tensor obeys the constraint 

^L = ^'Sta. (2.20) 



III. EVOLUTION AND CONSTRAINT EQUATIONS FOR THE BIANCHI II 
BRANE-WORLD 

A. Orthonormal Frame Formalism 

We shall adopt the orthonormal frame formalism as developed by Ellis and MacCallum 



p3| and further developed in ||2^. We shall assume that there is a group of motions 
acting simply transitively on three dimensional spacelike hypersurfaces. We will choose an 
orthonormal frame of vector fields {ca} and align eo to the fluid 4-velocity u, i.e., eo = u. 
Essentially we are assuming that the velocity field of the matter source is non-tilted with 
respect to the normal vector field of the hypersurfaces generated by the group of motions. 
Therefore, the triad of spacelike vectors {ea} at each point spans the tangent space of the 
group orbits and hence coordinates can be chosen such that the commutation functions are 
functions of t only. 

The commutation relation between eo and eo, yields the kinematic quantities associated 
with the vector field u. That is 

[eo,ea] = loaGb 



and jP^=-a/-H6^^-eJ(uj' + Q' 



where 

%a = Ua, and = -a^- - ^<^a - e^s 

Where H can be interpreted as the expansion of the fluid, 0"^/? can be interpreted as the 
rate of shear tensor of the fluid, iia can be interpreted as the acceleration of the fluid, 
and uJa is a measure of the vorticity of the fluid. The quantity Qa can be interpreted as 
the angular velocity of the spatial tetrad {ea} with respect to a Fermi-propagated spatial 
frame. However, with the assumptions made thus far both ita = and cUq, = 0. 



The vector fields generate a Lie Algebra with commutation functions 7^^ that can be 
decomposed in the following way 

[e„,e^] = -fipGs (3.1) 



{eap^n^'^ + aj^ - a^Sl^ (3.2) 



where riap and are functions of time t only. One can use a time dependent spatial 
rotation to diagonalize Wq,^ = diag(?T,i, ^2, ^3) p3|, |2^. The Jacobi identity [ei, [62,63]] = 
(and permutations) imply that Uapo,^ = 0, motivating the classification of a G3 group of 
motions into two classes, Class A where a" = 0, and Class B where a" 7^ 0. These two 
classes be be further sub-classified to yield the nine Bianchi types. 

In this paper we are interested in a set of Class A models and in particular the Bianchi 
type II models. For the Bianchi type II models only one of the najj is different from zero. 
Here we assume that rii = n, and 112 = = 0. In addition, we shall assume that the spatial 
triad {6^^} is Fermi-propagated which implies that f^Q, = 0. In choosing a Fermi-propagated 
triad, we have effectively determined the effective energy flux, that is, q^°^^^ = 0. This then 
imposes the condition that Qa = 0. In addition, in choosing a Fermi-propogated frame, one 
now has the freedom to rotate the spatial frame so that cTq/j = diag((7ii, (T22, ctss) ||23|. Since 
there is no way to determine Vab from observations on the brane, we impose the condition 
that our brane world model should be consistent with a Friedmann-LeMaitre cosmological 
brane model by assuming that Vai3 = 0. The equations describing the evolution of this class 



of models can be found in equations 1.90-1.100 of p4|. The Einstein Field equations and 
Jacobi identities yield 



H = -H^- -a^ - — (pt°t^i + 3p'"'^') 
3 6 


(3.3) 


2 2 

an = -3Han - -n 


(3.4) 


(J22 = -3Ha22 + 


(3.5) 


(T33 = -3i^'o-33 + 


(3.6) 


n = {-H + 2an)n 


(3.7) 


4 


(3.8) 


defined in equations p.l4-2.15| and we have introduced the constant 



where andp*°*'^^ . 

. Note, not all of the above equations are independent in particular an + (722 + (T33 = 0. 
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From the conservation equations on the brane, (see equations p.l8| and |2.19|) we obtain 



f> = ~3-yHp 

= -3H<j) - A;/t1/(0) 



(3.9) 
(3.10) 



Furthermore from the Bianchi identities on the brane (see equations |2.2U| and |2.11| ) we obtain 
(see 0), 

U = -AHU. (3.11) 



B. Dimensionless Formalism 

To facihtate the analysis of these models and to compare with previous work we now 
choose new variables of the form 



pf 



1 (722 + (^33 

2 H 



1 a22-q'33 

2V3 H - 



A 



3i/2' 



3 K(f) 

23H' 



N 



1 n 



(3.12) 



K 



and a new time variable 



6Aif2 2^ 
dt _ I 



V 



2^ 



The system of equations ( p.3| - p.ll]) becomes (where prime denotes differentiation with 
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respect to r). 



E'+ = J:+{q - 2) + (3.13) 

= S_(g-2) (3.14) 

N' = A^(g-4S+) (3.15) 

= Qpfi2q-i3-f-2)) (3.16) 

= 2^]A(g+ 1) (3.17) 

= l^A (2g - (37A - 2)) (3.18) 



^' = ^(g - 2) - ^A;$ (3.20) 
$' = 2$(g + 1 + ^A;^) (3.21) 



where 



q = 2J:l + 2J:l + ^^^Qpf~Q^ + ^^^^^Qx + ^u + 2^^-<^, 

where the generahzed Friedmann equation becomes 

(]t°tai = + + + + + $ = 1 - - - A^l (3.23) 

In addition, the equation for H decouples and becomes 

H' = -{l + q)H (3.24) 

We have now determined the equations describing the evolution of the Bianchi II brane- 
world model. The resulting equations are suitable for a qualitative analysis using techniques 
from dynamical systems theory. In general, if we let X = S_, A^, f^A, ^a, ^w, ^, 
then the system of equations ( |3.13| ), can be interpreted as X' = F(X) where F : X C — »• 
M^. Please see review of dynamical systems theory in and other texts such as p9|, ^ for 



a more complete description of dynamical systems techniques apphed to systems of ordinary 
differential equations. 
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IV. QUALITATIVE ANALYSIS OF THE BIANCHI II BRANE WORLD MODELS 



A. Invariant Sets 

We note that equation ( p.23[ ) is a constraint or conservation equation that we rewrite as 



G(X) = 1 - - - iV^ - - - - - $ - = 0. (4.1) 

Hence the phase space of the dynamical system S is an eight dimensional subset of M^, 
where 

5 = {XgM^|G(X) =0}. 

Here, one is able to use ( |4.1|) to eliminate the variable I^a globally from the dynamical system 
( p.l3| ), which results in a system of eight differential equations. 



The dynamical system ( 3.13|) is invariant under the transformation —N, therefore 



without loss of generahty we restrict the analysis to the set > 0. The set = {X G 
iS|A^ > 0} represents Bianchi H models, where as, the set = {X G S\N = 0} represents 
the Bianchi I models (if + ^ 0) or the zero curvature Friedman-Lemaitre models (if 
+ =0). 

The dynamical system ( p.l3|) is also invariant under the transformation S_ — S^, 
therefore without loss of generality we restrict the analysis to the set S_ > 0. 

The evolution equation for Qu implies that the surface Qu = ^ divides the phase space 
into three distinct regions, W+ = {X G 51^^^ > 0}, W° = {X G = 0} and = {X G 

s\nu < 0}. 

With the assumption of the weak energy condition (i.e., Ppf + + > 0), it can be 
shown that in the invariant sets and that 

< s^, s_, iv, npf, ^A, nx, nu, $, < i, 

that is, and W° are compact subsets of S. Unfortunately, the invariant set U~ is not 
compact in the variables we have chosen. One can however choose different dimensionless 
variables (|3.12|) to remedy this problem (see van den Hoogen and Abolghasem |28[| ). 
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There are various invariant sets associated with the matter content. We define six sets 



= {X G SlQpf = 0, $ = 0, ^ = 0}, 
W'^ = {X G = 0,$ = 0,^ ^ 0}, 

= {X G s\npf = 0, $ ^ 0, ^ ^ 0}, 
= {X G s\npf ^ 0, $ = 0, ^ = 0}, 

= {X G ^ 0,$ = 0,^ ^ 0}, 

+17+'^ = {X G 5|r]p/ ^ 0, $ 7^ 0, ^ 7^ 0}, 

where the notation is interpreted as 

(value of Qpf)^{value of <1>, value of NE") 

One additional invariant set of interest is the invariant set QTZ = {X G iSI^a = 0,^u = 
0}. One is interested in whether this set is the future attractor of all brane- world cosmologies. 

B. Monotonic Functions 

In the set U^, Q\ > 0, it can be shown that g > — 1. Therefore we easily see that 
monotonically increases from its lower limit = to its upper limit ^a = 1- Furthermore, 
with the existence of this monotonic function, we can rule out the existence of any closed or 
periodic orbits in the full eight- dimensional phase space, S. That is, any closed or periodic 
orbits (if they exist) must lie on the lower- dimensional boundaries of S. 

The seven-dimensional set U^, Q\ > lies on the boundary of S. If S_ ^ (or Qpf ^ 0) 
it can be shown that g > — 1 and again ^a monotonically increases. Again the existence of 
this monotonic function implies that there do not exist any closed or periodic orbits in this 
set. Note the only invariant subset of Qa > in which q = —1 (i.e, the only invariant 
subset of W° in which = 0) is S+ = S_ = = fipj = ^ = $ = and r^A + = 1- 

With the arguments made in the previous two paragraphs, the only sets in which a 
periodic and/or closed orbit might exist is the seven-dimensional set Q\ = and the five- 
dimensional set Qpf = Qu = S_ = 0,Q\ > 0. In all other cases we observe that the 
future asymptotic behaviour characterized by lim^^oo = 1 (this also implies through 
the constraint equation (|3.23|) that limT-_,ooX = [0,0,0,0,1,0,0,0,0]), the DeSitter Model 
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(see local analysis in next section), while the past asymptotic behaviour is characterized by 
limT-^_oo = 0, cosmological models with a four dimensional cosmological constant that 
is equivalently equal to zero. 

C. Equilibrium Points and Local Behaviour 

The equilibrium points can be classified in to one of the six invariant sets *■ ^^l^ ' \ depend- 
ing on the matter content. The order of the coordinates is S_, A^, ^Ipf, I^a, ^a, ^u, ^])- 
In what follows, for those points located in the Q\ = invariant set, the eigenvalue found 
in brackets { } is the eigenvalue associated with direction. 

1. Vacuum, 



Unfortunately the dynamical system (|3.13|) is not differentiable everywhere, and in par- 



ticular, it is not differentiable in a neighborhood of this invariant set. However, since the 
system is continuous, we are able to find the equilibrium points of the dynamical system. 
To determine the stability of each equilibrium point in this set, we change the variables 
{flpf, $) to spherical coordinates (p, 9, 0) and analyze the resulting system. This is done 
in Appendix 0. 

DS; De-Sitter [0, 0, 0, 0, 1, 0, 0, 0, 0]. 

m; Isotropic Braneworld Solution [0, 0, 0, 0, 0, 1, 0, 0, 0]. 

R^; Robertson Walker Dark Radiation Solution [0, 0, 0, 0, 0, 0, 1, 0, 0]. 

R^^] Bianchi II Dark Radiation Solution [i, 0, i, 0, 0, 0, |, 0, 0]. 

K] Kasner (Vacuum) Solution [cos(6'), sin(6'), 0, 0, 0, 0, 0, 0, 0] where —n < 9 < n. This 
point is part of the Kasner surface discussed below. 

2. Massless Scalar Field, ^Q^'^ 

/C; Kasner (Massless Scalar Field) Solution [sin((y9) cos(6'), sm{(p) sin(6'), 0, 0, 0, 0, 0, cos((y9), 0] 
where -n < 9 < n and -f < < f . Or + + = 1. At this equilibrium 
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point the deceleration parameter q — 2. The eigenvalues of the linearization restricted 
to the surface G'(X) = are (Qa eliminated) 

-6, 0, 0, 2, 2 - 4E+, 6 - 37, 6 + {6} 

The zero eigenvalues correspond to the fact that this is a two-dimensional set of equi- 
libria. We note that this equilibrium set has at least one negative and one positive 
eigenvalue, and hence this point is always a saddle, even in the JIa = invariant set. 
However, in the invariant set D,x — 0, a portion of this set becomes a local source. 

3. Massive Scalar Field, 

SF^; Robertson- Walker Scalar Field Solution (Power Law Inflationary Solution) 
[0, 0, 0, 0, 0, 0, 0, — 1 — ^] At this equilibrium point the deceleration parameter 
q — ^ — 1. The eigenvalues of the linearization restricted to the surface G'(X) = are 
(Qa eliminated) 

-k\ e -A,k'- 37, l(k' - 6), ^-(k' - 6), l(fc2 - 6), 1(A;^ - 2), {k'} 

We observe that this point has both negative and positive eigenvalues and that the 
point is not in the physical phase space for A;^ > 6. In the Qa = invariant set, 

• if /c^ < 2 and 1 < 7 < 2 then this point is stable. 

• if 2 < /c^ < 37 and 1 < 7 < 4/3 then this point has a 1 dimensional unstable 
manifold. 

• if 37 < A;^ < 4 and 1 < 7 < 4/3 then this point has a 2 dimensional unstable 
manifold. 

• if 4 < A;^ < 6 and 1 < 7 < 4/3 then this point has a 3 dimensional unstable 
manifold. 

• if 2 < /c^ < 4 and 4/3 < 7 < 2 then this point has a 1 dimensional unstable 
manifold. 

• if 4 < A;^ < 37 and 4/3 < 7 < 2 then this point has a 2 dimensional unstable 
manifold. 
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• if 37 < /c^ < 6 and 4/3 < 7 < 2 then this point has a 3 dimensional unstable 
manifold. 

We therefore conclude that for < 2, the power law inflationary solution is an 
attr actor in the JIa = invariant set. 



SF''; Bianchi II Scalar Field Solution [2^, 0, .J^^^g^, 0, 0, 0, 0-^, 36^^^]. 
At this equilibrium point the deceleration parameter q — 8^^^. This point only exists 
for parameter values in the range of 2 < A;^ < 8. The eigenvalues of the hnearization 
restricted to the surface G'(X) = are (Qa eliminated) 



-i Q c c k^-8 o 7fc^-32 



{k' - 8) ± ^{k' - 8)2 + 12(fc2 - 8)(fc2 _ 2) ,{18p^} 



We observe that this point has both negative and positive eigenvalues and therefore 
this point will always be a saddlepoint. In the JIa = invariant set 

• if 1 < 7 < 4/3 and 2 < k^ < 167/(6 - 7) then this point is stable. 

• ifl<7<4/3 and 167/(6 — 7) < k^ < 32/7 then this point has a one dimensional 
unstable manifold. 

• ifl<7<4/3 and 32/7 < k'^ then this point has a two dimensional unstable 
manifold. 

• if 4/3 < 7 < 2 and 2 < /c^ < 32/7 then this point is stable. 

• if 4/3 < 7 < 2 and 32/7 < k'^ < I67/ (6 — 7) then this point has a one dimensional 
unstable manifold. 

• if4/3<7<2 and 167/(6 — 7) < k'^ then this point has a two dimensional 
unstable manifold. 

RSF^; Robertson Walker Radiation-Scalar Field Scaling Solution 
[0, 0, 0, 0, 0, 0, — gp-]. At this equilibrium point the deceleration pa- 
rameter q — 1. The eigenvalues of the hnearization restricted to the surface G'(X) = 
are (Qa eliminated) 

-4, -1, -1, 4 - 37, 1, -^{k ± VA:^ + 16(4-F)), {4} 
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This point exists in the set C/+UC/° only if k'^ > 4. We observe that this point has both 
negative and positive eigenvalues and therefore this point will always be a saddlepoint. 
In the JIa = invariant set 

• ifl<7<4/3 then this point has a 2 dimensional unstable manifold. 

• if4/3<7<2 then this point has a 1 dimensional unstable manifold. 

RSF^^; Bianchi II Radiation-Scalar Field Scahng Solution [|, 0, |, 0, 0, 0, gp] 
At this equilibrium point the deceleration parameter q = I. The eigenvalues of the 
linearization restricted to the surface G(X) = are eliminated) 

-4, -1,4 -37, {4} 

(^{2k^) ± y(i2)^^^^2^ (23A;2 - 64) ± ^{23k^ - 64)2 _ 64A;2(7A;2 - 32) ^ 

This point exists in the set U only if k"^ > 32/7. This point is a saddle in the 
full phase space. However, in the invariant set fl\ = 0, 

• ifl<7<4/3 then this point has a one dimensional unstable manifold. 

• if4/3<7<2 then this point is stable. 



4. Perfect Fluid Models, 

F°; Robertson Walker Perfect Fluid Solution [0, 0, 0, 1, 0, 0, 0, 0, 0] At this equilibrium point 
the deceleration parameter q = |(37 — 2). The eigenvalues of the hnearization re- 
stricted to the surface G(X) = are {fl\ eliminated) 

-37, ^(7 - 2), ^(7 - 2), ^(7 - 2), 37 - 4, ^(37 - 2), 37, {37} 

We observe that this point has both negative and positive eigenvalues and therefore, 
this point will always be a saddlepoint. However, in the invariant set = 0, 

• ifl<7<4/3 then this point has a 2 dimensional unstable manifold. 

• if4/3<7<2 then this point has a 3 dimensional unstable manifold. 

F^^; Bianchi II Perfect Fluid Solution [|(37 - 2), 0, 1^3(37 - 2)(2 - 7), ^(6 - 
7), 0, 0, 0, 0, 0] At this equilibrium point the deceleration parameter q — |(37 — 2). 
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The eigenvalues of the hnearization restricted to the surface G(X) = are (Qa ehmi- 
nated) 

-37, 37 - 4, 1(7 - 2), 1(7 - 2), 37, {37}, 
-I (2(2 - 7) ± V[2(2-7)F + 2(2-7)(7-6)(37-2)) 

We observe that this point has both negative and positive eigenvalues and therefore, 
this point will always be a saddlepoint. However, in the invariant set D,\ — 0, 

• ifl<7<4/3 then this point has a 1 dimensional unstable manifold. 

• if4/3<7<2 then this point has a 2 dimensional unstable manifold. 

5. Perfect Fluid with a Massless Scalar Field, 

There are no equilibrium points in the interior of this invariant set. 

6. Perfect Fluid with a Massive Scalar Field, "'"ri"'"'^ 

MSF^; Robertson Walker Matter-Scalar Field Scahng Solution [0,0,0,1 - 
||, 0, 0, 0, — ^^^^^^] At this equilibrium point the deceleration parameter 
q — i(37 — 2). The eigenvalues of the linearization restricted to the surface G'(X) = 
are (Qa eliminated) 

-37, 1(7 - 2), 1(7 - 2), 1(37 - 2), 37 - 4, {37}, 
|((7-2)±^(7-2)2 + 87(7-2)(l-fj)) 

This point only exists in the physical phase space when k"^ > 37. We observe that this 
point has both negative and positive eigenvalues and therefore, this point will always 
be a saddlepoint. However, in the invariant set JIa = 0, 

• ifl<7<4/3 then this point has a 1 dimensional unstable manifold. 

• if4/3<7<2 then this point has a 2 dimensional unstable manifold. 

MSF^^; Bianchi II Matter-Scalar Field Scaling Solution [|(37 - 
2),0,|V3(37-2)(2-7),l - + 1^(2 - 37), 0, 0, 0, At this equi- 
librium point the deceleration parameter q — |(37 — 2). The eigenvalues of the 
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linearization restricted to the surface G{X.) — are {Q\ eliminated) 

-37, 1(7 - 2), 37 - 4, {37}, -I (^A ± ^A^ + ± VB^+C)^ 

where 
A = 2(2 - 7) 
S = 247(1^-1) + 3(2 -7)2 
C = 647(37 - 2)(f!:+7-6) 

This point only exists in the physical phase space when 1 — p + ^(2 — 37) > 0. We 
observe that this point has both negative and positive eigenvalues and therefore, this 
point will always be a saddlepoint. However, in the invariant set Q,\ — 0, 

• ifl<7<4/3 then this point is stable. 

• if4/3<7<2 then this point has a 1 dimensional unstable manifold. 

V. OBSERVATIONS 

The early time behaviour of the brane- world cosmologies under consideration is no longer 
characterized by a Kasner type singularity. It is found that the initial singularity is repre- 
sented by the isotropic brane-world solution. It is also observed that this initial singularity 
is very much dominated by the kinetic energy of the scalar field. It is worthy to note that 
the Kasner solutions are not stable to the past in the brane-world scenario, however, they 
continue to have an unstable manifold having a large dimension, and are indeed stable to 
the past for General Relativistic cosmologies. 

It is no surprise that the De Sitter solution represents the late time behaviour of any cos- 
mology with a positive four dimensional cosmological constant A. The early time behaviour 
of any cosmology with a positive four dimensional cosmological constant is asymptotic to 
a regime with an equivalent four dimensional cosmological constant equal to zero. That is, 
the early time behaviour of cosmologies having a positive four-dimensional cosmological con- 
stant fl\ > 0, can be determined by investigating the behaviour of the governing dynamical 
system restricted to the Qa = invariant set. 

We observe that the late time behaviour in the = invariant set is dependent upon 
the curvature of the models under consideration. For example, the late time behaviour of 
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the Bianchi II cosmological models is one of four possible behaviours/solutions. If A;^ < 2 
then the isotropic power-law inflationary solution is the stable attractor. If A;^ > 2 then the 
late time behaviour is one in which the models fail to isotropize, and have either a scalar 
field source, a combination of a scalar field and an ordinary matter source, or a combination 
of a scalar field and a dark radiation source depending upon the values of k and 7. (See table 
|V|). On the other hand, the late time behaviour of the Bianchi I cosmological models is also 
one of four possible solutions. If < 2 then the isotropic power-law inflationary solution 
is the stable attractor. If A;^ > 2 then the late time behaviour is one in which the models 
continue to isotropize, and have either a scalar field source, a combination of a scalar-field 
and ordinary matter source, or a combination of a scalar field and a dark radiation source 
depending upon the values of k and 7. (See table for summary details). 

For the most part, it becomes very obvious that the existence of the scalar field in these 
brane-world cosmologies affects the initial dynamics when compared to the perfect fluid 
models studied in |^ . The scalar field also plays an important role in the future dynamics 
of those models with a zero cosmological constant, that is, the scalar field is a dominant 
feature in the future dynamics of both the Bianchi type II models and the Bianchi type 
I models. Two solutions of particular interest, the radiation-scalar field scaling solutions 
{RSF^^ , RSF^) are shown to be stable for particular values of 7 and k. Since neither of 
these solutions exist in general relativistic cosmologies, and since both of these cosmologies 
are due to bulk effects acting on the brane, we have the conclusion, that in general (contrary 
to popular belief), these brane-world scenario's are not necessarily asymptotic to general 
relativistic cosmologies to the future. 

The qualitative analysis for the bifurcation values 7=1 and 7 = 4/3, and for the various 
values of k has not been completed. However, this analysis is not expected to yield any new 
observations. Additional numerical calculations do not reveal anything unexpected. See 
figures for phase portraits. Note, the analysis presented here also does not address the 
case in which the local energy density due to non-local effects from the free gravitational 
field is negative, Qu < 0- If < 0, then these traditionally ever-expanding models could 
potentially recollapse due to the negativity of this dark radiation. This analysis is presently 



under consideration in [ES 
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TABLE I: Sinks and Sources for the Bianchi type II models: Set represents the invariant set 
containing the basin of attraction/repulsion, and the Dimension is the dimension of the basin of 
attraction 



Equihbrium Point 


Sink/Source 


Set 


Dimension 


Conditions 


m 


Source 




8 


7 > 1 


DS 


Sink 




8 




JC 


Source 


nx = o,N+ 


7 


E+ < 1/2,^ < -V6/k 


SF^ 


Sink 


= 0, 


7 


e<2 


SF^^ 


Sink 


Qa = 0, N+ 


7 


1 < 7 < 4/3, 2<P < 167/(6 - 7) and 










4/3 < 7 < 2,2 < < 32/7 


MSF^^ 


Sink 


^A = 0, N+ 


7 


l<7<4/3,fc2 > 167/(6-7) 


RSF^^ 


Sink 


^A = 0, 


7 


4/3 < 7 < 2, P > 32/7 




Sink 


^A = 0, N+ 


7 


4/3 < 7 < 2 


RSF'^ 


Sink 


Qa = 0, 


6 


4/3 < 7 < 2, > 4 


SF° 


Sink 


Qa = 0, 


6 


1 < 7 < 4/3, < 37 and 










4/3 < 7 < 2,/c2 < 4 


MSF^ 


Sink 


^A = 0, N° 


6 


1 < 7 < 4/3, fc2 > 3^ 




Sink 


flA = 0, 


6 


4/3 < 7 < 2 
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Appendix: Equilibrium Points in the Vacuum Invariant Set 



We define spherical coordinates as follows, 



tan^e = -^7 (5.1) 

2 7 ^pf 

COS d) — 



The evolution equations for (p, 6, (p) are 

p' = p(q+l _ ^cos^^-Ssin^^cos^^ ) , (5.2) 



9' = sm9{ ScosO + ^kpsm(l) ] , (5.3) 
sin cos 0(^ — 3 cos^ 9 ) (5.4) 



where 



g = 2S^ + 2S2 + ^^y^p2 cos^ - + fi^^ + gin^ 0(2 cos^ ^ - sin^ 9) + ^^^^ '^ Qx 
and 

7a = 27 cos^ 0^ + 4 sin^ cos^ 9. 
The constraint equation (|4.1| ) becomes 

+ nj, + nx + ^u + ^l + ^- + N^ -^^- = (5.5) 

and we have < 9 < n and < < 7r/2 since both $ > and Qpf > 0. 

If we define Y = (p, ^, 0, S+, S_, A^, i^A, ^w) then we obtain a new dynamical system 
Y' = F(Y). Where the dynamical equations for S+, S_, A^, Qx, are the same as those in 
equation ( p.l3|) with the change of variables (^.1|) , and is given by solving equation ( ^.5|) 



for ^A- Note, we indicate the eigenvalue corresponding to the ^a direction with brackets 
{}• 

The geometrical interpretation of the analysis of this new system requires some care, 
as each equilibrium point will have a number of different coordinatizations in this new 
coordinate system. The dynamics in the rectangular coordinate system {Qpf, \E', $) can be 
obtained by determining the dynamics along invariant rays or directions (the intersections of 
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invariant planes {9' = 0) with invariant cones = 0)) in the spherical coordinate system. 
We quickly observe here that 6 = 7t/2, (p = tt/2 does not correspond to an invariant direction 
of the dynamical system Y' = F(Y). 

DS; De-Sitter 

The equilibrium point corresponding to the De Sitter solution in the original variables 
is represented by six different equilibrium points in the spherical coordinates defined 
above: DS = [0, 6', 0, 0, 0, 0, 0, 0] where 6 = 0,7r/2,7r and = 0,7r/2. The eigenval- 
ues — 3, — 3, — f , — 4, have an eigenspace spanned by T,^,T,_, N,Qii. The eigenvalues 
associated with the matter are Xp = —37/2, A^^ = —67 when = and Xp = —3, 
Xq^ = —12 when = 7r/2 except at the point (0 = 7r/2,0 = 7r/2) where the eigenvalues 
are Xp = 0, Aq^ = 0. The eigenvalues associated with the two angular coordinates are: 
Ag- = 3 when = O,7r,0 = 0,7r/2) and A^ = -3(2 -7)/2 when = 0, vr, = 0), and 
X^ = 3(2-7)/2 when = 0, vr, = 7r/2); Ag = -3 when {6 = 7r/2,0 = 0,7r/2) and 
A_^ = 37/2 when {6 = 7r/2,0 = 0), and A^ = -37/2 when {6 = 7r/2,0 = 7r/2). The 
behaviour of the angular coordinates 0, 6 is illustrated in Figure We observe that 
all the eigenvalues not associated with the angular coordinates are negative, there- 
fore, this point is a local attractor. The dynamics in the (f) — 9 plane indicate that 
— * 7r/2 and ^^^7r/2asr^ 00. This implies that this equilibrium point is strongly 
attracting along this direction, that is, this equilibrium point strongly attracts along 
the VLpf = \E' = direction in the original variables. This local analysis, together with 
the global analysis in section |1VB| , implies that the De Sitter solution is a global at- 



tractor for a wide class of Bianchi II, Bianchi I, and zero curvature Robertson Walker 
cosmological models with a cosmological constant. 

m; Isotropic Braneworld Solution 

The equilibrium point m in the original variables is represented by six different equi- 
librium points in the new variables, m = [0,^^,0,0,0,0,1,0] where 6 = 0,7r/2,7r 
and = 0,7r/2. For 9 = 0,7r/2,7r and = the eigenvalues associated with the 
non-angular coordinates are {67}, 37/2, 3(7 — 1),3(7 — 1),37 — 1,2(37 ~ 2). For 
6 = 0,71 and = 7r/2 the eigenvalues associated with the non-angular coordinates 
are {12}, 3, 3, 3, 5, 8. For 9 = n/2 and = tt/2 the eigenvalues associated with the 
non-angular coordinates are — 3, — 3, — 1, 0, {0}, — 4. The eigenvalues of the angular 
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coordinates are the same as in De-Sitter case above (see Figure |l|). We observe that 
this point is always a source for 7 > 1 and strongly repel away from the Qpf = $ = 
direction. As models evolve backwards in time, it is the kinetic energy of the scalar 
field that is the dominant component of the total energy density. For 7 = 1, numerical 
analysis reveals that this point remains a source (see Figure 0). 

R^; Robertson Walker Dark Radiation Solution 

The equilibrium point in the original variables is represented by six different equi- 
librium points in the new variables, R^ = [0, 6', 0, 0, 0, 0, 0, 1] where 6 = 0,71/2, it 
and (p = 0, 7r/2. The eigenvalues — 1, — 1, 1, {4} have an eigenspace spanned by 
T,^,T,_, N,Qii. The eigenvalues associated with the matter are Ap = (4 — 37)/2, 
Xq^ =4 — 67 when = and Ap = —1, A^^^ = —8 when = 7r/2 except at the point 
{9 = 7r/2,0 = 7r/2) where the eigenvalues are Ap = 2, Aq^ = 4. The eigenvalues of 
the angular coordinates are the same as in De-Sitter case above. We observe that the 
eigenvalues not associated with the angular coordinates are both positive and negative, 
therefore, this point is a saddle point. If 7 < 4/3 then the point has a 3 dimensional 
unstable manifold in the rectangular coordinate system. (Note: If 7 < 4/3, then 
Ap = (4 — 37) /2 > along the invariant direction = but Ap = — 1 < along the 
invariant directions ^ = 0,0 = 7r/2 and ^ = 0,0 = vr/2.) If 7 > 4/3, then the point 
has a 2 dimensional unstable manifold. 

R^^; Bianchi II Dark Radiation Solution 

The equilibrium point R^^ in the original variables is represented by six different 
equilibrium points in the new variables, R^^ = [0, 6^, 0, |, 0, |, 0, |]. The eigenvalues 
-l,-l/2(l±V6i),{4}have an eigenspace spanned by S+, S_, A^, ^lu- The eigenvalues 
associated with the matter are Ap = (4 — 37)/2, A^^ = 4 — 67 when = and Ap = — 1, 
Xq^ = —8 when = 7r/2 except at the point {6 = 7r/2, = tt/2) where the eigenvalues 
are Ap = 2, A^^ = 4. The eigenvalues of the angular coordinates are the same as in De- 
Sitter case above. If 7 < 4/3 then the equilibrium point has a 2 dimensional unstable 
manifold in the rectangular coordinate system. If 7 > 4/3, then the equilibrium point 
has a one dimensional unstable manifold. We also observe that if 7 > 4/3 then this 
point is a stable attractor in the Q\ = invariant set. 
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FIG. 1: Behaviour of the equilibrium points in the (p — 9 plane. 
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FIG. 2: Qualitative behaviour of trajectories when 7 = 1, A; = 1 and k = 2. 

Phase Portrait in the projected 1'-4> plane. (y= 1 and k = 1) Phase Portrait in the projected "P-^) plane. (y= 1 and k = 2) 
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FIG. 3: Qualitative behaviour of trajectories when 7 = 7/6 < 4/3, k = 1 and k = 2. 



Phase Portrait in the projected T-* plane. (y= 7/6 and k = 1) Phase Portrait in the projected 'i'-* plane. (y= 7/6 and k = 2) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



27 




28 




29 



FIG. 6: Qualitative behaviour of trajectories when 7 = 5/3 > 4/3, k = 3. 



Phase Portrait in the projected 'i'-* plane, (y = 5/3 and l< = 3) 
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